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The super symmetric approach in the form of second order intertwining re- 
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lations is used to prove the exact solvability of two-dimensional Schrodinger 
equation with generalized two-dimensional Morse potential for oq = —1/2. 
This two-parametric model is not amenable to conventional separation of 
variables, but it is completely integrable: the symmetry operator of fourth 
order in momenta exists. All bound state energies are found explicitly, and 
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all corresponding wave functions are built analytically. By means of shape 
invariance property, the result is extended to the hierarchy of Morse models 
with arbitrary integer and half-integer values a\. = — (k + l)/2. 
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^ ■ Introduction. - Very limited list of exactly solvable one-dimensional models is known 

up to now: harmonic oscillator, singular harmonic oscillator, Coulomb, Morse, Poschl- Teller 
potentials, and some others [1]. As for complete solution of the Schrodinger equation in 
two-dimensional space, no regular approach is known. In practice, only the conventional 
separation of variables [2] provides the sole method to solve two-dimensional stationary 
quantum problems [3]. One of the fundamental problems of modern Quantum Mechanics 
is to enlarge as much as possible the variety of exactly solvable multi-dimensional (d > 2) 
systems. 

The supersymmetric approach has provided a powerful impulse for new development in 
analytical studies both in one- dimensional [6] and multi-dimensional [7] Quantum Mechanics. 
But to date, even among two-dimensional potentials, not amenable to standard separation 
of variables [8], very few results on partially solvable models were obtained. By means of 
supersymmetry, several two-dimensional models were found [9] - [12], for which a part of their 
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spectra and wave functions were found analytically (partial solvability). It is necessary to 
stress that all constructed systems are completely integrable, since the symmetry operators 
of fourth order in momenta were built for them. Two specific supersymmetric techniques 
- the special SUSY-separation of variables and two-dimensional shape invariance - allowed 
to investigate two-dimensional singular generalizations of Morse [9] and Poschl-Teller [11] 
potentials. The first method explores the most important constituent of SUSY algebra - 
SUSY intertwining relations. The second method is the natural generalization of the one- 
dimensional shape invariance approach [13], which provided an elegant algebraic algorithm 
for solution of a class of Schrodinger equations. In one-dimension, the class of shape invariant 
potentials coincides practically with the known variety of exactly solvable models [1]. 

Among all possible two-dimensional intertwining relations with supercharges of second 
order in derivatives, a subclass exists, where one of intertwined Hamiltonians is amenable to 
standard separation of variables due to specific choice of parameters of the model. Such situ- 
ation was studied already for the complex two-dimensional Morse potential [15], which is free 
of any singularities. In the present Letter, we investigate the model, which is much more in- 
teresting physically - the real singular two-dimensional Morse potential with parameter value 
a = —1/2. Namely, we prove that, being completely integrable, this two-parametric model 
is exactly solvable as well, i.e. all its eigenvalues and eigenfunctions are known analytically. 
After that, we use two-dimensional shape invariance of the model to enlarge the class of 
exactly solvable two-dimensional systems to a hierarchy of Morse potentials with parameter 
values ak = —(k + l)/2; k — 1, 2, 



Morse potential with a = —1/2. - Two-dimensional generalization of Morse potential 
[9, 10] is defined [14] as: 



where a is an arbitrary real parameter, and U (xj) are one-dimensional Morse potentials: 



V(x) 




1) 



+ U(xi) + U(x 2 ), x = (xi, x 2 ), x± = xi± x 2 , 



(1) 



U( Xi ) = A(e 



—2ctXi 



<le~ ax % i = 1, 2 A = Const > 0; a = Const > 0, 



(2) 



The Hamiltonian H(x) 



A ( 2 ) + V(x) with potential (1) and the partner Hamiltonian 
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H{x) with potential: 

tW . a 2 a(2a + 1) rW , . . . . 

^ = si nh 2 (^) + ^ + (3) 

are intertwined [8] - [10] 

H(x)Q + = Q + H{x)- Q-H(x) = H{x)Q- (4) 

by the second order supercharges: 

Q ± = 4<9 + <9_ ± 4aa<9_ ± 4aa coth ^^<9 + + 4a 2 a 2 coth - 

2 2 
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One can notice that V and V differ from each other by a coefficient in front of singular 
term only, this property being the origin of shape invariance of the model [16]. By construc- 
tion, each intertwined Hamiltonian obeys the symmetry operator of fourth order in momenta 
[9, 10], which is expressed in terms of supercharges Q ± . Specifically, the operator R = Q + Q~ 
commutes with Hamiltonian H(x), and operator R = Q~Q + - with H{x). 

For the particular choice of parameter a = —1/2, the partner Hamiltonian (3) is sim- 
plified essentially - it allows the conventional separation of variables. Moreover, each of 
one-dimensional problems, obtained after separation, is exactly solvable one-dimensional 
Morse potential (2). The discrete spectrum of this model is well known [17]: 
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e n = -a s n , s n = n-->0; n = 0,1,2,..., (6) 

a 2 

and wave functions are expressed in terms of degenerate hypergeometric functions [18]: 

rjn(xi) = exp(-^)(&) s "F(-n, 2s n + 1; & = exp(-ax i ). (7) 

2 a 

Due to separation of variables, the quantum problem with the Hamiltonian H(x) from (3) 
is exactly solvable. Its energy eigenvalues are: 

E nt m E m ^ n 6 n + 6 m , (8) 

being two-fold degenerate for n ^ m. The corresponding eigenfunctions can be chosen as 
symmetric or (for n 7^ m) antisymmetric combinations: 



Our aim here is to solve completely the quantum problem for H{x) with potential (1) for 
a = —1/2. The main tool are the SUSY intertwining relations (4), providing [10] connections 
between spectra and wave functions of partner Hamiltonians, which are almost isospectral 
(up to zero modes and singular properties of Q^). 

In general, we may expect three kinds of levels of H(x). 

(i) . The levels, which coincide with (8). Their wave functions can be obtained from (9) 
by means of intertwining relations (4). 

(ii) . The levels, which were absent in the spectrum of H(x), if some wave functions of 
H(x) are simultaneously the zero modes of the supercharge operator Q~ . Then the second 
intertwining relation in (4) would not give any partner state among bound states of H{x). 

(iii) . The levels, which were also absent in the spectrum of H{x), if some wave functions 
of H(x) become nonnormalizable functions after action of operators Q~ . 

We will analyze these three classes of possible bound states of H(x) one after another. 

(i). The first supersymmetric intertwining relation in (4) gives immediately the two-fold 
degenerate wave functions of H(x) with energies (8), as ^> E nm = Q + ^E nm - Using the explicit 
formulas (5) for the supercharge Q + and for one-dimensional Schrodinger equation (2), one 
can rewrite the symmetric (antisymmetric) wave functions in the form: 

= Q + *k m = - + D*k, m , (ii) 

where the differential operator D is defined as: 

D = T^-Ei + 6 + 266(<% + %)]. (12) 
Since this operator is singular for £i = £ 2 , normalizability of ^ E nm depends crucially on 

A S 

the behavior of ^E nm on the line £i = £2- The symmetric wave functions on this line are 
*fL, m (fi»6 = fi) ~ Vni&rimiZi), and the Y have no zero multipliers - £ 2 ), which could 
compensate the singularity of operator D. Therefore, the wave functions ^E nm = Q + ^E nm 
are certainly nonnormalizable, and they do not correspond to the physical bound states. 

Vice versa, action of the operator D in (10) on the antisymmetric function l ea ds 
to nonsingular result. To observe this property, one has to use the explicit form (7) of one- 
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dimensional eigenfunctions. After separation of exponential multiplier exp[— (£1 + ^)/2], 
the rest of function ty E is a polynomial in £1,^2- Due to antisymmetry, this polynomial 
vanishes for £1 = £2, since it contains the multiplier (£1 — £2), compensating the singularity 
in D. Moreover, careful study shows that owing to the interplay between two terms in the 
r.h.s. of (10), the symmetric wave functions ^f nm ~ (£1 — at the singular point £1 = £ 2 . 

It is possible to investigate efficiently the normalizability of wave functions ^f. by the 
indirect algebraic method as well. One can check straightforwardly that for a = — 1/2 the 
operator R = Q~Q + can be rewritten as: 

R = (/ii(xi) - h 2 {x 2 )^ + 2a 2 [h l {x 1 ) + h 2 (x 2 )j + a\ (13) 

where hi{xj) = — A i + U(xi) are one-dimensional Schrodinger operators with Morse potential 
(2). Acting by this operator on the antisymmetric function , one obtains: 

R*i n Jx) = [(e n -e m ) 2 + 2a 2 E ntm + a 4 ]^i n Jx) = 

= A(n ~ m? ~ l][(a„ + s m ) 2 - lM n Jx) = r n ^ En Jx). (14) 

This relation helps to find the norm of wave functions : 

W*kJ 2 = «. m I Q~Q + I *1,J = rn, m \\*tj 2 , (15) 

where we have used that Q ± are mutually conjugate on the considered space. We stress that 
off-integral terms in (15) at the singular point x 1 = x 2 disappear due to vanishing behavior of 
normalizable wave functions *ff s (x), mentioned above (see below Eqs.(17),(18) and nearby). 

By definition, ^ E m in (15) are zero identically for m — n. In its turn, the wave functions 
^E nn±1 are annihilated [15] by the action of Q + , since the corresponding r n m in (14) vanish 
for m = n ± 1. It follows from (14), that the norms of all other ty En m (with | n — m \ > 2) are 
finite and positive. Therefore, the first class of levels of H(x) consists of nondegenerate energy 
levels (8) with normalizable symmetric wave functions ^J; from (10) with | n — m \ > 2. 

(ii). This class of possible bound states of the Hamiltonian H(x) with a = —1/2 con- 
sists of the normalizable zero modes of the supercharge Q~ from (5). In [9] the variety of 
normalizable zero modes of the Hermitian conjugated supercharge Q + was studied in detail 
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for arbitrary values of a. In particular, it was shown that these zero modes are normalizable, 
and no fall to the center occurs for the specific range for values of a. According to expressions 
(5), Q~ can be obtained from Q + by replacing a — > —a. Therefore, completely analogous 
analysis of zero modes of Q~ will lead to the analogous, but positive, interval for a : 



A 1 
a e (- + — — , +00), 



(16) 



^4 ' 4^/2 

which, for certain, does not contain the value a = —1/2. This means that no normalizable 
bound states of this class exist for H. 

(iii). We have to study an opportunity that operator Q~ destroys normalizability of 
some eigenf unctions of H. It could occur due to singular character of Q~ at x± — x 2 - It is 
convenient to choose the polar coordinates £, f in the vicinity of this line in the plane (£1, £ 2 ) : 

£i=£cosy?; £ 2 =£sinv?; ip = tt/4 - 9; 9 ~ 0. 



In these variables the potential V(x) reads: 
V = a 2 



cos 29 1 9 . 3. r-^ 

+ 7^ o)v^cos0 



sin^ 



Thus, the Hamiltonian if (a:) effectively acts as: 

H = -A (2) ~ -a 

and two kinds of behavior of its eigenf unctions are possible: 



2 de 9* 



9 2 



m ~9 2 or ~ _1 . 

Only the first one is normalizable at # ~ with the measure dx\dx 2 = d^d9 /£ cos 20. 
In the same vicinity, the supercharge operator Q~ acts as: 



(17) 



(18) 



Q-~« 2 (^-l)(<% + i). 



This operator is not able to transform normalizable wave function of (18) to nonnormalizable. 
Therefore, the third class of possible wave functions \I/ of H(x) does not exist too. 

Summing up the above analysis, we found that the spectrum of two-dimensional Hamil- 
tonian H(x) with a = —1/2, which is not amenable to separation of variables, consists only 
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of the bound states with energies (8) for \n — m\ > 1. This spectrum is bounded from above 
by the condition of positivity of s n , s rn in (6): n,m < y/A/a — 1/2. The corresponding wave 
functions are obtained analytically, according to (7), (9), and (10). 

Hierarchy of Morse potentials with a k = —(k + l)/2. - In what follows we denote 
the exactly solvable Hamiltonians H{x), H(x) investigated above, as H(x;a ), H(x;a ). 
The motivation is that we will consider the Hamiltonians H(x; a k ), H(x; at) with arbitrary 

negative integer and half-integer values a = a k = —(k + l)/2; k — 0, 1, 2, We will prove 

that all these Hamiltonians are also exactly solvable. 

The important property is that due to identity a k -i(2a k ^i — 1) = a k (2a k + 1) : 



This means that the following chain (hierarchy) of Hamiltonians can be built: 
H(x; a ) -T- H(x; a ) = H(x; a ± ) -j- H(x; a ± ) = . . . + H(x; a k -i) = H(x; a k ) 4- H(x; a k ), (20) 

where the sign -f- between two Hamiltonians denotes their intertwining by Q ± {ai). 

Since the Hamiltonian H(x; a ) (and therefore H(x; ai)) was shown to be exactly solvable, 
we can use the intertwining relations from (20) to prove exact solvability of H(x; a±) too. 
Considering this procedure for the general length of the chain (20), we will illustrate results 
for the first section k — 1 of the chain, with the Hamiltonian under study H(x; a x = —1). 

It was shown above, that bound states of the Hamiltonian H(x;a ) = H(x;ai) are 
described by the symmetric functions ^% mn {x\ oq) = ^% mn (x\ a\) of (10) and energy eigen- 
values (8) with \n — m\ > 1. Therefore, the functions 



*l im (f ; ai) = Q + (ai)*| njm (f; ai) = Q + ( ai )Q + (a )^i n Jx; a ); \n - m\ > 1 (21) 



(if they are normalizable) are the wave functions of the partner Hamiltonian H(x; ai). The 
only exclusion concerns possible zero modes of Q + (a 1 ). It was shown in [9], that the wave 
functions of H(x; a k ) may coincide with zero modes of supercharge Q + (a k ) for energy levels: 



H(x; a k -i) = H(x; a k ); k = 1,2,.... 



(19) 



2 / 2 \ 2 2 

-En,n+fc+i = ~2a {2a k s n + s n ) - Aa a k 



2/2 i 2 

a (s n + s n 



n+k+l 
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In particular for k — 1, operator g + (ai), acting on the wave functions ^ E n ,„± 2 (^i a i) °f 
ai), may annihilate them, not admitting to the spectrum of H(x; ai). 
Analogously in the general case, up to zero modes of Q + the functions 

^E n , m (x; a k ) = Q + (a k )^ En ^(x; a k ) = Q + (a k )Q + (a k ^) . . . Q + (a )^i n Jx; a ) (22) 

(again, if normalizable) are the wave functions of H(x; a k ) with energies E n>rn = — a 2 (sf l +s^ n ). 
The symmetries of wave functions (22) alternate and depend on the length of chain (20). 

In order to control the situation with possible zero modes of Q + and normalizability 
of the states (22), we will use the following algebraic trick. For general value of a k , the 
identity between Hamiltonians H(x; a k ) = H(x; a k -\) leads to identity (up to a function of 
the Hamiltonian itself) between their symmetry operators R(a k ) and R(a k ^i). This relation 
can be derived straightforwardly by algebraic manipulations with result: 



2H(x; a fe _i) + a 2 (2k 2 + 2k + 1) 



(23) 



Q~(a k )Q + (a k ) = Q + (a k ^)Q~(a k ^) + a\2k + 1) 

In the case k = 1, Eq.(23) can be used to transform the norm of the function (21): 

ll^, m (^«i)H 2 = (^l im (f;a )|g-(a )g-(a 1 )g + (a 1 )g + (ao)|^l im (f;a )) =r„, m (24) 
(r n , m + 6a 2 E n<m + 15a 4 )||^ m (a )|| 2 = a 4 r n>m [(m - n) 2 - 4][(s n + s m ) 2 - 4]||^ m (a )|| 2 , 

where definitions (14) for r„ im were explored. From Eq.(24), the (real) functions ^ n ,n±2(x] a±) 
vanish identically for m = n ± 2, being zero modes of Q + (a\). Therefore, energies £ , „, n +2 
are absent in the spectrum of H(x; a x ). Vise versa, due to the same Eq.(24) the norms of all 
other *& n>m (x; ai) with \n — m\ > 2 are positive and finite. 

This analysis can be generalized to arbitrary integer values of k. The proof is performed 
step by step, by using Eq.(23) for indices k, k — 1, . . . , 1, in calculation of the norm: 

11*^(3 «fc)f = (*i, m (3 ao)\Q-(ao) ■ ■ ■ Q-(a k )Q + (a k ) . . . Q + (a )|*i >m (f ; a )>. (25) 
After some algebra, one obtains that 

11*^(3 «fc)H 2 = <*i, m (3 a )\R(a ) (Rfa) + r ) (r(oo) + T + T^j . . . 

. . . (Rfa) + r + . . . + r fc _x) \^i n Jx; oo)); (26) 
T ; = a 2 (2l + l)[2H(a ) + a 2 {2l 2 + 21 + 1)]. 



Inside the matrix element (26), R(ao) can be replaced by its eigenvalues (14), and sums 
(r„ + r\ + . . . + Vi) by: 

-2(4 + slW + 2) + (. + l) 4 -l". 
Thus, the norm (25) can be factorized: 

«*) II 2 = a 4 *!!*^^; ao) - "») 2 - (* + 1) 2 ][(*» + *™) 2 - + 1) 2 ]- (27) 

i=i 

One can conclude that the norm vanishes for \n — m\ < (k + 1), and it is finite positive 
for other n,m. This means that the spectrum of H(x; a k ) includes all energy levels E n>m of 
H(x; ao) with \n — m\ > (k + 1) only. All others disappear due to zero modes of Q + (cii). 

Similarly to the analysis (ii) above, no normalizable wave functions of H(x; a k ) can be 
annihilated by the operator Q~(a k ), since the values a k = —{k + l)/2 lie again outside of 
the interval (16), permitted for their normalizability. The class (iii) of possible bound states 
of H{x ; a k ) is also empty. To prove this fact, one has to consider the behavior of operators 
H, Q~ and normalizable wave functions ^ at the point of singularity 9 ~ : 

Q-{a k ) ~ a 2 fa - (k + 1)) (d e + *±1). 

Similar to the case k — 0, operators Q~(a k ) can not destroy normalizability of ^ E n , m {%] a k)- 
Summarizing the obtained results, the spectra of Hamiltonians H(x; a k ) are not de- 
generate. They consist of the bound states with energies E njfn , given by (8) with indices 
\n — m\ > k + 2. Their wave functions ^E nm {x; a k ) are given by Eq.(22) in terms of degen- 
erate (confluent) hypergeometric functions. Like in the case ao = —1/2, the discrete spectra 
are bounded from above by the conditions (6): n,m < \[Ala - 1/2. 
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